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Abstract. Electromagnetic waves and fluids have locally conserved mechanical 
properties associated with them and we may expect these to exist for matter 
waves. We present a semiclassical description of the continuity equations relating 
to these conserved properties of matter- waves and derive a general expression for their 
respective fluxes. 



It has long been known that electromagnetic waves have associated mechanical 
properties. Maxwell described the radiation pressure exerted by light pQ and Poynting 
identified the momentum carried by light [2] and that angular momentum was associated 
with circular polarisation [Hj. The energy, momentum and angular momentum of 
electromagnetic waves are particularly important as these quantities are all conserved. 
Associated with the density of each of these quantities is a flux density (or current 
density) which is related to the density via continuity equations [I]. 

The creation of the first Bose-Einstein condensates (BECs) in dilute atomic gases 
in 1995 E] brought about the possibility of creating coherent matter waves. A BEC 
in a trap is analogous to an electromagnetic field in a laser cavity, in that macroscopic 
occupation of a single quantum state occurs. Just as a laser cavity is a source of coherent 
light, a BEC can be a source of coherent matter waves. To continue this analogy we 
assert that matter waves will have associated mechanical properties and present a general 
theory of the local mechanical properties of matter waves. We give, as examples, the 
densities and flux densities of the matter wave particle number and energy. The same 
theory applies to their momentum and angular momentum. 

The local conservation of particles is a well established concept in fluid mechanics. 
It states that in any volume element of a fluid, no matter how small, the rate of change 
of particles in an element of fluid is equal to the rate at which the particles flow into the 
surface of the element. This principle is made quantitative in the following equation [Hj; 

»/ v/ riV = -/^.dS, (1) 

where p is the particle density, v is the velocity of the fluid and dS is a surface element 
with direction normal to the surface enclosing the volume V. If we define the particle 
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current density (or flux density) J = pv, then we can easily obtain the continuity 
equation for particles in a fluid 

| + V.J = 0. (2) 

Similar continuity equations hold for other conserved quantities [S] such as the energy 
density 



dw „ 

sT + v 




0, (3) 



where e is the internal energy per unit mass and P is the pressure of the fluid, and the 
momentum density 

^ + A. (PSij+pv , V]) = . (4) 

Here the summation convention has been introduced, so the indices (i and j) take the 
values 1, 2, 3 and a repeated index is understood to be summed. This convention will 
be assumed for the remainder of the paper. 

In a case in which a globally conserved quantity does not obey the continuity 
equation of the form (J2J), there will be a source of that quantity. An example of this is a 
fluid under the influence of an external force, which is a source of momentum. In such a 
system it is necessary to modify the equation of continuity to allow for the source term. 
We would expect that the rate of change of momentum in an element of volume would 
be equal to the rate at which momentum flows into the surface plus the rate at which 
momentum is created in the element. This reasoning leads us to a modified continuity 
equation of the form 

where n ij? is the {i,j} component of the momentum flux density tensor and Si is the 
source density for the i component of momentum. 

We now turn our attention to the semiclassical description of matter waves. 
Consider an observable A having a corresponding Hermitian operator A. The density 
A(R) of A at a point R will obey the generalised local continuity equation 

^A{R) + ^% {A \R)=& A \R), (6) 

where Xi are the components of R, % is the flux density of the observable and is 
the source density term. We define the Hermitian operator associated with the density 
of A at R to be 

A{R) = \{A,8{t-R)}, (7) 

where {, } denotes the anticommutator and S(r — R) is the Dirac delta function. From 
the Heisenberg equation of motion we can show that the density will evolve according 
to 
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We see that the first term on the right hand side will vanish if A is conserved. It is 
therefore reasonable to associate this term with the source term in (JSJ and the last 
term with the flux density. In order to proceed in the calculation we use the standard 
Hamiltonian for a particle in a potential V{r) 



H = —ViVi + V(r) 
2m 



where V,- 



-Jj^-. On evaluating the last term in 




we find that 

1 f d ~ 
= 5 L^(r-R) 



(9) 



(10) 



which takes the form of a conservation equation (jHJ). This allows us to identify the flux 
density operator 



% iA) (R) 




(11) 



We can associate the operator {pi/m,8(r — R)} with the velocity density and hence 
see that T^ A \R) is a suitably symmetrised product of the velocity density and the 
observable. This operator is the flow of the density of the observable, that is its flux 
density. The forms of the flux density and source density for the observable A will come 
from the expectation value of (JTUJ). In the Heisenberg picture, the matter wave will be 
described by a wavefunction ip(r). The expectation value of the conservation equation 
for this state is 



d 



^A(R) + V R ^ A \R)=S^(R 



(12) 



which is the continuity equation for the expectation values of the relevant operators. 
Here we have introduced the notation A(R) = (A(R)). The flux density expectation 
value is 

T {A) < 



R) 



4m 



d 3 r^*(r)^ (A) (i?)^(r) 



r=R 



and the source density expectation value is 



S (A) (R) = ±{-([H,A]l>) m 1> + p([H, A] ^) } 



r=R 



(13) 



(14) 



As an example we consider two quantities which obey continuity equations, specifically 
the particle number and energy. Following Landau 0, the operator corresponding to 
particle density in a semiclassical theory is represented by the Dirac delta function, so 
in (J7J) we put the operator A equal to the identity operator 1. From (fT3*j) the flux density 
of particles at position r can be calculated as 



% {P) (R) 



-Im {rVM\ r =R 



(15) 



where Im indicates the imaginary part. This is the form of the well-known quantum 
mechanical probability flux Jj ^U]- The source term vanishes as the identity operator 
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commutes with the Hamiltonian. This is a reasonable result as there is no particle 
source present in the system. Energy in quantum mechanics is represented by the 
Hamiltonian operator, so we can calculate the energy flux density from (jlHj) . If we use 
the Hamiltonian © we obtain the energy flux density at position R, which is 



The source term vanishes as the Hamiltonian commutes with itself. 

In this letter, we have derived general flux density and source density operators for 
conserved quantities in matter waves within the semiclassical description. This theory 
will be important when considering the deposition of these quantities during interactions 
of matter waves with other objects. An example of this is the detection of matter waves 
themselves, which requires one to consider the particle flux density [TT] . One problem is 
that the flux density operator, expression (jllj) . might contain an arbitrary divergenceless 
operator in addition to the terms present. It seems likely, as with Poynting's theorem, 
that this freedom will have no physical consequences j!2j . Also of interest is the form of 
the second-quantised version of the theory presented here. We will address these issues 
elsewhere. 
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4m 2 



Im {t/A7 2 V^* + (ViV*) V 2 V>} 



r=R 



+ V(R)% W (R).(16) 



